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Abstract 



O 

^SJ ' The total gluon helicity in a polarized proton, measurable in high-energy scattering, is shown 

i-t • to be the large momentum limit of a gauge-invariant but non-local, frame- dependent gluon spin 



E X A± in QCD. This opens a door for a non-perturbative calculation of this quantity in lattice 
QCD, and also justifies using free-field expressions in the light-cone gauge as physical observables. 



The total gluon helicity AG* = J dxAg{x) in a longitudinally polarized proton is an 
important physical quantity that has motivated much experimental effort to measure in 
high-energy scattering \^]^. It helps to understand, among others, how the helicity of a 
fast-traveling proton is composed of partons' helicity and orbital angular momentum |5|- 
|9|. The factorization theorems in quantum chromodynamics (QCD) indicate that AG is a 



matrix element of a complicated light-cone correlation operator of the gluon fields [10[, and 
has a siniple physical interpretation only in the light-cone gauge A~^ = natural for parton 
physics |5|. This state of understanding makes computing AG in lattice QCD infeasible, 
and raises the fundamental question about the gauge invariance of the gauge particle spin 



in a bound state 11 



In this paper, we report a breakthrough in understanding the physics of AG, and 
correspondingly leading to a practical way to its calculation in lattice QCD. We find that 
AG can be obtained by boosting a matrix element of the gluon spin operator — E x A_j_ — to 
the infinite momentum frame (IMF), where A± is the transverse part of the gauge potential. 
This operator was first proposed in Ref. 1^ as the gauge-invariant gluon spin, but has been 



criticized as physically uninteresting because of its frame dependence [13[. The physics 



behind the IMF limit we propose here goes back to the well-known Weissacker- Williams 



equivalent photon picture for high-energy scattering [ll|. But of course, there are subtleties 
in taking the IMF limit. In particular, the matrix element has a singular dependence on 
the bound state momentum in perturbation theory as it approaches infinity. Moreover, the 
anomalous dimension oi E x Aj_ does not coincide with that of the non-local light-cone 
correlation in the factorization theorems [ij]. Therefore, we will provide a well-defined 
procedure, or a matching condition, for the limiting procedure. In particular, we will show 
in a one- loop example how to obtain AG from the IMF limit of a frame-dependent, time- 
independent matrix element. This example helps to demonstrate that a non-perturbative 
AG can be recovered from a frame- dependent lattice matrix element of an Euclidean space 
operator E x A_i_. 

The difficulty in understanding and calculating AG is easy to appreciate. Through QCD 
factorization, it has been shown that AG is a matrix element of a non-local operator involving 



light-cone correlation |10| . 

AG = j dx^ j '^e--''^^- {PS\Fr{n^''\r,^)K:,mPS) , (1) 

where \PS) is a proton plane- wave state with momentum P^ and polarization S*^, F"^ = 

(1/2) e'^f^'^'Ff,^, and £(^,0) = Pexp[-ig jf ^+(r/-,Ox) dr]-] with A+ = T''A+ is a light- 
cone gauge link defined in the adjoint representation. The light-front coordinates are defined 
as ^^ = (^° ± ^^)/\/2. It is usually difficult to see the above operator as the gluon spin 
or helicity. However, in the light-cone gauge, A'^ = 0, the whole operator collapses into 
E X A — the textbook definition of the gauge particle spin ll| — which is known to be gauge 
dependent, and AG can be regarded as the number of gluon partons with helicity in the 
direction of the proton helicity minus that with opposite helicity. Because of the explicit 
presence of the real time in ^~, one cannot evaluate the above expression in lattice QCD. 
An early attempt to get the gluon helicity on lattice was to calculate the matrix element of 



F^yF^^ [15[, but there is no demonstrated connection between this and AG. 

To find the physics of AG without committing to the light-cone gauge, let us examine 
the operator structure a bit further. For simplicity, we first consider the U(l) gauge theory 



(quantum electrodynamics, or QED) so that the gauge hnk is absent. Carrying out the 
integration over the longitudinal momentum, the gauge-invariant photon "spin" operator 
becomes: 



ff--(0) = ^ /— /^ IdCd'Ue-'^''''^^'-'^-^^^^ 



xj (27r)3 
dk'^d'^kx 
(27r)3 



^xP-'A\r,^±)-^klA^{d±) ^.+(0) 



A\k^,k^)--^A^{k+,k±) F,+(0) 



k+ 



^(O)x(A(O)-— VA+(r 



(2) 



where V"*" = d/dC,~ , -E* = -F*"*", k~^ = xP~^, and the third component of a vector is interpreted 
in the usual sense of the cross product. The ^~ coordinate in A~^ is taken to after operating 
with the inverse derivative. 

The above operator is just the IMF limit of E' x A±, where E^ = F*°, and A^j_ = {6^^ — 
V^V-'/V^)/!-^ is the transverse part of the gauge potential and is invariant under gauge 
transformation. The rule of taking the IMF limit of an operator is as follows: For any 
vector V^^ define V^ = (V° ± V^)/^/2. If the boost is along the 3-direction, then the 
components of the vector go like V~^ — )■ ^"'"A, V" — )■ V~ /A, and Vj_ — )■ Vj_. Thus V^ — )■ 
(V"'")^A^ and V ■ A — )■ V'^A'^A'^ for the leading components. Using these rules, one finds 
that E X A± ^ E X {A- (1/V+)VA+), which is exactly the above operator in Eq. ([2]). 

It is known that in electromagnetic theory the vector potential can be uniquely separated 
into the longitudinal and transverse parts, A = Ay + A±, and the transverse part is gauge 
invariant [16|: Given A, A± can be uniquely constructed as a functional of A with an 
appropriate boundary condition. Thus E x A±_ is a gauge-invariant operator, and can be 
regarded as the gauge-invariant part of the gauge particle spin 12|. E x A± is a non-local 
operator in that it depends on the gauge potential over all space. 

It is important to realize that separating A and E into longitudinal and transverse parts 
is in general not a physically meaningful thing to do. In the first place, the physics of 
E is to apply force to electric charge and there is no charge that responds separately to 
E\\ and E^. Second, in a different frame, one sees different transverse and longitudinal 
separations, and therefore the notion has no Lorentz covariance [13|. As we shall see, the 
frame- dependence of both parts is dynamical, and cannot be calculated without solving the 
theory. However, there are two exceptions where the separation is meaningful. The first 
case concerns the radiation field [17|. For free radiation, by separating out the unphysical 
degrees of freedom, one simplifies the quantization procedure significantly. In the laser 
beam, this separation allows one to talk about the gauge-invariant photon spin and orbital 



angular momentum [18|. The second case is the IMF, which is our interest here. In the 
IMF, E'li -C E_j_, and the electromagnetic field can be regarded as free radiation. This has 
been recognized long ago by Weissacker and Williams, in the name of equivalent photon 
approximation [ll|. Only then, E x A± can be understood as a physical quantity, where E 
can also be replaced by E±. 

Therefore, the photon helicity measurable in high-energy scattering can be calculated as 
the IMF limit of a matrix element of the static operator E x A±. To calculate the matrix 
element of the time-independent, albeit non-local, operator is a standard practice in lattice 
QCD. It will be dynamically dependent on the momentum of the external particle [13[. In 
fact, the dependence is singular in the leading order perturbation theory: In the IMF limit. 



the matrix element diverges. To obtain the physically interesting finite light-cone matrix 
element, one has to find a matching condition, which we will come to after establishing a 
similar connection in QCD. 

The case for QCD is a bit more complicated. Separating A into longitudinal and trans- 
verse parts requires generalizing the observations in QED to similar conditions in QCD, 
which has been considered long ago [19| (see also Ref. 12|). Clearly, we would like to have 
A± transform covariantly under gauge transformation: 



A^ -^ U{x)A^U^{ 



x) 



(3) 



where A±_ = TaA"^, so it is easy to construct gauge-invariant quantities with A±. Second, 
we require A\\ to produce null magnetic field, as it does in QED. This condition is 19] : 



aM|" - d'A]:" - gr'^A^M 







(4) 



which is a nonlinear equation to solve for A* as a functional of A^. Moreover, the transverse 

(5) 



part of the gauge potential satisfies a generalized Coulomb condition [19 

aM^ = ig[A\A^j_] . 



We can then go through a similar derivation as in the QED case and find: 



(6) 



where the inverse derivative acts on everything after it, and takes the ^~ coordinate in the 
gauge link to 0. 

It is a bit involved to show that the expression in the parenthesis above is indeed A± in 
the IMF. Since A± = A — A\\, we just need to solve for A\\. After solving Eq. (^ and Eq. ([5]) 
order by order in g, we find that Aj^ vanishes, and thus At = A'^ in the IMF. Substituting 
this into Eq. @, we obtain a first-order inhomogeneous linear equation for Al: 



a+Aj" - gf^^A'^^^'A^f = d'A-^^" . 
Its solution is easy to construct as a geometric series expansion: 



Al 



V^ 



1 + 



'ZgA^ 



1 



-igA" 



1 



ab 



(9M^ 



v+y V V 

By commuting c^M"*" systematically to the front of the expression, one finds: 



(7) 



(8) 



(9) 



where the coordinate ^'~ in the gauge link is taken to ^~ after operating with the inverse 
derivative. 

Alternatively, one can multiply a gauge link C on both sides of Eq. ([7]), and find after 
some manipulations: 

d+iA'fC'"^) = (9M+'")£"'^ . (10) 



It is then straightforward to see that the solution for An'" is formally given by Eq. (Q. 

Clearly A||(0) is just the part subtracted from A°-{0) in the gluon helicity operator 
in Eq. ([6]). Therefore, we established the same conclusion in QCD that the gluon helic- 
ity operator is the IMF limit of the gauge-invariant, non-local gluon spin operator {E x A±)^. 

Now we show that the matrix element of S*™^ = {E x Aj_)^ depends on the choice of 
frames dynamically. By "dynamically" we mean that the frame dependence cannot be 
obtained from Lorentz transformation, and is a function of dynamic details. Let us consider 
the example of photon spin in a free electron state \p, s). A simple perturbative calculation 
of Fig. [1] yields. 




FIG. 1: Matrix element oi E x Aj_ in an asymptotic electron state. 
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P 
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ra^ 



u{p, s)T,^u{p, s) . (11) 



where p'^ = (p°,0,0,p^), D = 1/e — 7£;-|-ln47r-|-ln(/i^/?7i^), and m is the mass of the electron 
used to regularize the infrared (IR) divergence. The result has a non-trivial dependence on 
the electron momentum p^, which makes its physical interpretation less straightforward. In 
particular, the result is divergent in the IMF limit. The ultra-violet (UV) part of the matrix 
element is consistent with that found in Ref. [ij]. Our result also shows that the statement 
in Ref. [20| about the frame independence is incorrect. 

However, the measurable photon helicity can be obtained by the same matrix element 
by going to the IMF limit first before the loop-momentum integration. In this limit, the 
external momentum dependence is dropped out, and the matrix element becomes: 



{p,s\S^^\p,s) 



Ur. 



An 



i3D + 7) u{p,s)J:^u{p,s) . 



(12) 



This is exactly the same as that computed using the factorization expression in Eq. ([T]), or 



using Ex A (notice the full A), in the light-cone gauge A~^ = [2]|. The UV property of the 
matrix element is the same as that derived by Altarelli and Parisi (AP) [22] • In fact, in the 
original derivation of AP evolution equation, a finite frame result was used to obtain parton 
physics in the IMF limit. Our result also indicates the claim that gluons carry only about 
1/5 of the nucleon momentum in Ref. 12|] is incorrect, because the matrix elements of the 
quark and gluon momentum operators were calculated in the finite momentum frame. We 
have verified that if they are boosted to the IMF, we can get the standard mixing matrix 
in Ref. |23|. Meanwhile, it should be pointed out that the result in Ref. 
because they used the light-cone gauge condition. 



24 is standard 



The above calculation shows that the IMF and UV limits are not exchangeable. However, 
since the IR part is the same, the difference is a perturbatively calculable quantity. This 
turns out to be the key for a non-perturbative computation of the gluon helicity. Since 
lattice QCD cannot handle the real time dependence, a direct calculation of Eq. ([T]) is 
infeasible. However, one can get the same matrix element by studying the matrix element 
of S"™^ = E X A_j_ as a function of external momentum P. The largest momentum attainable 
on a lattice is of order 1/a, with a being the lattice spacing. On the other hand, the matrix 
element also has UV dependence on 1/a, and can be calculated in perturbation theory. Thus 
we can calculate the matrix element at the largest momentum ~ 1/a, matching the results 
of the two different limits in a perturbative way. For instance, in the above one-loop example 
with a finite e, one can match the two results by setting: 

16 

\np^ = {D + \nm^) + 21n2. (13) 

o 

The result will become more accurate as the renormalization scale is asymptotically large. 
On the lattice, D + Inrri^ is replaced by — Ina^, and the matching condition becomes 
In(pa)^ = const. We will explore the issue of the lattice calculation in more detail in a 
separate publication |25| . 

At last, it is useful to consider the frame dependence of the angular momentum sum 



rule which has recently been strongly advocated in the literature [12|, |20|, |26| . The primary 
goal is to find a simple free-field form of the angular momentum decomposition so that 
the individual parts have simple physical interpretation. The closest gauge-invariant form 
involves the expression |l2l ]: 

J = (fx tp'^-I^i) + / d^x tp'^x X -(V - igA\\)ilj 

+ [ d^x EaXAl+ [ d^x Eixx VA'f . (14) 



This result is frame dependent and not physically interesting in general. However, in the 
IMF one has A^f{^~) = (l/V+)(V^A+'^£'"'(^'",^~)); the above decomposition is the same 
as the Jaffe-Manohar result in the light-cone gauge A~^ = |5|. Therefore, this serves 
to justifying that the light-cone gauge is the natural choice in the IMF, where free-field 
expressions such a.s E x A attain physical significance. In particular, the corresponding 
matrix elements are physically measurable. 

To conclude, we have shown that the total gluon helicity measured in high-energy scat- 
tering is the IMF limit of a matrix element of a gauge-invariant operator. This limit does 
not commute with the UV limit in quantum field theory and therefore the two operators 
have different anomalous dimensions at first sight. However, they can be related through 
a matching condition. This allows an otherwise infeasible lattice QCD calculation of light- 
cone correlations. We have also explained why free-field theory expressions in the light-cone 
gauge are physically meaningful, as they correspond to the IMF limit of gauge- invariant, 
but non-local expressions in interacting theories. 

We thank J. -W. Chen for useful discussions about the gauge-invariant gluon spin op- 
erator. This work was partially supported by the U. S. Department of Energy via grants 

6 



DE-FG02-93ER-40762 and a grant (No. 11DZ2260700) from the Office of Science and Tech- 
nology in Shanghai Municipal Government, and by the National Science Foundation of China 

(No. 11175114). 



[1] The COMPASS Collaboration, Phys. Lett. B 633, 25 (2006); Phys. Lett. B 647, 8 (2007); 
arXiv: [hep-ex] 0802.3023 (2008); Phys. Lett. B 676, 31 (2009); Phys. Lett. B 718, 922 (2013); 
Phys. Rev. D 87, 052018 (2013). 
[2] The HERMES Collaboration, JHEP 1008, 130 (2010). 

[3] The PHENIX Collaboration, Phys. Rev. Lett. 93, 202002 (2004); Phys. Rev. D 76, 051106 

(2007); Phys. Rev. D 79, 012003 (2009); Phys. Rev. Lett. 103, 012003 (2009); Phys. Rev. D 

83, 032001 (2011); Phys. Rev. D 84, 012006 (2011); Phys. Rev. D 86, 092006 (2012); Phys. 

Rev. D 87, 012011 (2013). 

[4] The STAR Collaboration, Phys. Rev. Lett. 97, 252001 (2006); Phys. Rev. Lett. 100, 232003 

(2008); Phys. Rev. D 80, 111108 (2009). 
[5] R. L. Jaffe and A. V. Manohar, Nucl. Phys. B 337, 509 (1990). 
[6] S. Brodsky, M. Burkardt, and I. Schmidt, Nucl. Phys. B 441, 197 (1995). 
[7] R. L. Jaffe, Phys. Lett. B 365, 359 (1996). 

[8] X. Ji, J. Tian, and P. Hoodbhoy, Phys. Rev. Lett. 76, 740 (1996). 
[9] P. Chen and X. Ji, Phys. Lett. B 660, 193 (2008). 
[10] A. V. Manohar, Phys. Rev. Lett. 66, 289 (1991). 

[11] J. D. Jackson, Classical Electrodynamics, (3rd Edition, John Wiley & Sons Inc. 1999). 
[12] X. -S. Chen, X. -F. Lii, W. -M. Sun, F. Wang, and T. Goldman, Phys. Rev. Lett. 100, 232002 

(2008); Phys. Rev. Lett. 103, 062001 (2009). 
[13] X. Ji, Phys. Rev. Lett. 104, 039101 (2010); Phys. Rev. Lett. 106, 259101 (2011). 
[14] X. -S. Chen, W. -M. Sun, F. Wang, and T. Goldman, Phys. Lett. B 700, 21 (2011). 
[15] J. E. Mandula, Phys. Rev. Lett. 65, 1403 (1990). 

[16] F. J. Belinfante, Phys. Rev. 128, 2832 (1961); F. Rohrlich and F. Strocchi, Phys. Rev. 139, 
B476 (1965); V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Quantum Electrody- 
namics, (2nd edition, Pergamon, Oxford, 1982); C. Cohen- Tannoudji, J. Dupont-Roc, and 
G. Grynberg, Photons and Atoms, (John Wiley & Sons Inc. 1989). 
[17] X. Ji, Y. Xu, and Y. Zhao, JHEP 1208, 082 (2012). 
[18] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, Phys. Rev. A 45, 8185 

(1992). 
[19] R. P. Treat, J. Math. Phys. 13, 1704 (1972). 
[20] M. Wakamatsu, Phys. Rev. D 83, 014012 (2011). 
[21] P. Hoodbhoy X. Ji, and W. Lu, Phys. Rev. D 59, 074010 (1999). 
[22] G. Aharelh and G. Parisi, Nucl. Phys. B 126, 298 (1977). 
[23] H. D. Politzer, Phys. Rep. 14, 129 (1974); D. J. Gross and F. Wilczek, Phys. Rev. D 9, 980 

(1974); H. Georgi and H. D. Politzer, ibid. 9, 416 (1974). 
[24] M. Wakamatsu, Phys. Rev. D 85, 114039 (2012). 
[25] X. Ji, J. -H. Zhang, and Y. Zhao, in preparation. 
[26] M. Wakamatsu, Phys. Rev. D 81, 114010 (2010). 



